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This paper is devoted to the explicit construction of optimal de- 
signs for discrimination between two polynomial regression models of 
degree n — 2 and n. In a fundamental paper, Atkinson and Fedorov 
[Biometrika 62 (1975a) 57-70] proposed the T-optimality criterion 
for this purpose. Recently, Atkinson [MODA 9, Advances in Model- 
Oriented Design and Analysis (2010) 9-16] determined T-optimal de- 
signs for polynomials up to degree 6 numerically and based on these 
results he conjectured that the support points of the optimal design 
are cosines of the angles that divide half of the circle into equal parts 
if the coefficient of x™" 1 in the polynomial of larger degree vanishes. 
In the present paper we give a strong justification of the conjecture 
and determine all T-optimal designs explicitly for any degree n G N. 
In particular, we show that there exists a one-dimensional class of 
T-optimal designs. Moreover, we also present a generalization to the 
case when the ratio between the coefficients of x n ~ and x" is smaller 
than a certain critical value. Because of the complexity of the opti- 
mization problem, T-optimal designs have only been determined nu- 
merically so far, and this paper provides the first explicit solution of 
the T-optimal design problem since its introduction by Atkinson and 
Fedorov [Biometrika 62 (1975a) 57-70]. Finally, for the remaining 
cases (where the ratio of coefficients is larger than the critical value) , 
we propose a numerical procedure to calculate the T-optimal designs. 
The results are also illustrated in an example. 
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1. Introduction. The problem of identifying an appropriate model in 
a class of competing regression models is of fundamental importance in re- 
gression analysis, and it occurs often in real experimental studies. It is widely 
accepted nowadays that good experimental designs can improve the perfor- 
mance of discrimination, and several authors have addressed the problem of 
constructing optimal designs for this purpose; see Hunter and Reiner (1965), 
Stigler (1971), Atkinson and Fedorov (1975a, 1975b), Hill (1978), Fedorov 
(1980), Denisov, Fedorov and Khabarov (1981), Studden (1982), Fedorov 
and Khabarov (1986), Spruill (1990), Dette (1994, 1995), Dette and Haller 
(1998), Song and Wong (1999), Ucihski and Bogacka (2005), Wiens (2009, 
2010) among many others. In a fundamental paper, Atkinson and Fedorov 
(1975a) introduced the T-optimality criterion for discriminating between 
two competing regression models. As an example, these authors constructed 
T-optimal designs for a constant and a quadratic model. Since its introduc- 
tion, the problem of determining T-optimal designs has been considered by 
numerous authors; see Atkinson and Fedorov (1975b), Ucihski and Bogacka 
(2005), Wiens (2009), Tommasi and Lopez-Fidalgo (2010), among others. In 
order to demonstrate the benefits of the T-optimal design, we display, in Ta- 
ble 1, the simulated power of the T-test for the hypothesis Hq : #2,2 = #2,3 = 
in the cubic regression model rj(x,6) = S|=o^ 2 i x ' 7 on ^ ne interval [—1,1] 
(with standard normal distributed errors), where observations are taken ac- 
cording to two designs. The first design is the commonly used equidistant 
design with 12 observations at the four points, —1, —1/3, 1/3 and 1, respec- 
tively, while the second design is a T-optimal design, as considered in this 
paper, with 8 observations at the two points, —1,1 and 16 observations at 
the two points —1/2,1/2, respectively. We observe clear advantages (with 
respect to the power of the T-test) for the T-optimal design. 

Since its introduction T-optimal designs have found numerous applica- 
tions including such important fields as chemistry of pharmacokinetics; see 
Atkinson, Bogacka and Bogacki (1998), Asprey and Macchietto (2000), Ucihs- 
ki and Bogacka (2005) or Foo and Duffull (2011) among others. The T- 
optimal design problem is essentially a minimax problem, and, except for 
very simple models, the corresponding optimal designs are not easy to find 

Table 1 

Simulated power of the F-test in a cubic regression model 
'^2^ = o^2jX-' for the hypothesis of linear regression model for 
various values of #2,3 and different designs (82.2 = 0) 



02,3 





0.5 


1.0 


1.5 


2.0 


T-optimal 


0.051 


0.104 


0.301 


0.641 


0.896 


Equidistant 


0.053 


0.092 


0.218 


0.438 


0.638 
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and have to be determined numerically. In a recent paper, Dette and Titoff 
(2009) discussed the T-optimal design problem from a general point of view 
and related it to a nonlinear problem in approximation theory. As an illus- 
tration, designs for discriminating between a linear model and a cubic model 
without quadratic term were presented, and it was shown that T-optimal 
designs are, in general, not unique. 

Atkinson (2010) considered a similar problem of this type and studied 
the problem of discriminating between two competing polynomial regression 
models which differ in the degree by two. This author determined T-optimal 
designs for polynomials up to degree 6 numerically where the coefficient 
of x* 1 " 1 in the polynomial of larger degree (say n) vanishes. Based on these 
results he conjectured that the support points of the T-optimal design are 
cosines of angles dividing a half of circle into equal parts. 

The present paper has two purposes. In particular, we prove the conjec- 
ture raised in Atkinson (2010) and derive explicit solutions of the T-optimal 
design problem for discriminating between polynomial regression models of 
degree n — 2 and n for any n € N. Moreover, we also determine the T- 
optimal designs analytically in the case when the ratio of the coefficients 
of the terms and x n is sufficiently small. The situation considered in 

Atkinson (2010) corresponds to the case where this ratio vanishes, and in this 
case we show that there exists a one-dimensional class of T-optimal designs. 
To the best of our knowledge these results provide the first explicit solution 
of the T-optimal design problem in a nontrivial situation. Our results pro- 
vide further insight into the complicated structure of the T-optimal design 
problem. Finally, in the case where the coefficient exceeds the critical value, 
we suggest a procedure to determine the T-optimal design numerically. 

2. The T-optimal design problem revisited. Consider the classical re- 
gression model 



where the explanatory variable x varies in the design space X , and obser- 
vations at different locations, say x and x' , are assumed to be uncorrelated 
with the same variance. In (2.1) the quantity e denotes a random variable 
with mean and variance a 2 , and r/ is a function, which is called regres- 
sion function in the literature. We assume that the experimenter has two 
parametric models for this function in mind, that is, 



and the first goal of the experiment is to discriminate between these two 
models. In (2.2) the quantities 9\ and O2 denote unknown parameters which 
vary in compact parameter spaces, say 0i C M. mi and 02 C M m2 , and have 



(2.1) 




(2.2) 





4 



H. DETTE, V. B. MELAS AND P. SHPILEV 



to be estimated from the data. In order to find "good" designs for discrim- 
inating between the models rji and 772 , we consider approximate designs in 
the sense of Kiefer (1974), which are defined as probability measures on the 
design space X with finite support. The support points of an (approximate) 
design £ give the locations where observations are taken, while the weights 
give the corresponding relative proportions of total observations to be taken 
at these points. If the design £ has masses oj{, > at the different points X{ 
(i = 1, . . . , k), and N observations can be made by the experimenter, the 
quantities uJiN are rounded to integers, say rij, satisfying Yli=i n i = N, and 
the experimenter takes rij observations at each location X{ (i = 1, . . . , k). 

To determine a good design for discriminating between the models 771 
and 7)2 [Atkinson and Fedorov (1975a)] proposed in a fundamental paper to 
fix one model, say r\\ (more precisely its corresponding parameter 9\), and 
to determine the design which maximizes the minimal deviation between 
the model 7/1 and the class of models defined by 772 , that is, 

£*=argmax / (7/1(2;, 0i) - 772(2:, 9%)) 2 £(dx), 
where the parameter 9^ minimizes the expression 

6>2 = argmin / (771(2:, 0\) - 772(2;, 9 2 )) 2 ^(dx). 
6*2662 J x 

Note that 9\ is not an estimate, but it corresponds to the best approxima- 
tion of the "given" model 771 (-, 9\) by models of the form {772 (•, 9 2 ) \ 9 2 € 62} 
with respect to a weighted L2-norm. Since its introduction the T-optimal 
design problem has found considerable interest in the literature, and we refer 
the interested reader to the work of Ucihski and Bogacka (2005) or Dette 
and Titoff (2009), among others. In general, the determination of T-optimal 
designs is a very difficult problem, and explicit solutions are — to our best 
knowledge — not available except for very simple models with a few param- 
eters. In this paper we present analytical results for T-optimal designs, if 
the interest is in the discrimination between two polynomial models which 
differ in the degree by two. To be precise, we consider the case where the 
regression functions 7/1(2;, 9\) and r] 2 (x,9 2 ) are given by 

(2.3) 7/1(2:, 9 X ) = 0i,o + 0i,i2: + • • • + 0i,„_ 2 x n ^ 2 + 0i, n -i2: re - 1 + 9 lin x n 
and 

(2.4) 7/2(2:, 9 2 ) = 2 ,o + 02,12: + • • • + 9 2in - 2 X n ~ 2 , 

respectively, and the design space is given by X = [—1, 1]. In model (2.3) the 
parameter 0i is given by 0i = (0i,o 5 01,1, • • • > 0l,n-2, ^0l,ri) 0l,n) ■> where the 
ratio of the coefficients corresponding to the highest powers b = 9\ n —\/0\ n 
and the parameter 9\ n specify the deviation from a polynomial of degree 
n-2. 
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In the following discussion, we define 
fj(x,a,b,e ln ) = ?7i(a:,0i) - mix, 82) 

(2.5) 

= a + aixH ha n - 2 x n 2 + 8^ n {bx n 1 +x n ), 

where we use the notation cti = On — Q<i,% (i = 0, ... ,71 — 2); then the problem 
of finding the T-optimal design for the models 771 and m can be reduced to 

£* = argmax / {al + a\x^ h a*„ 2 £™~ 2 + #1 „(&x n_1 + x n )) 2 £(dx), 

where a* = (a\, . . . ,0^-2)^ ^ s a vector minimizing the expression 

a* = argmin / (fj(x,a,b,9i tn )) 2 ^(dx). 
a Jx 

It is now easy to see that for a fixed value of b = &i n -i/8in, the T-optimal 
design does not depend on the parameter 9± n . In the next section we give 
the complete solution of the T-optimal design problem if the absolute value 
of the parameter b = 0i, n — i/0in less or equal to some critical value. 

3. T-optimal designs for small values of |6| = |#i,n-i/#i,n|- Through- 
out this section we assume that the parameter b satisfies 



7T \ \ / / _ I TT \ \ 2' 71 " 



(3.1) |6| = |0i,,i_i/0i,„|<nl 1-cosl - 1 1 / I 1+cosl - 1 )=ntan^l — 
then it is easy to see that all points 

(3.2) tm =J 1+ m cos (™)J^ i=1 ,..., n , 

are located in the interval [—1, 1]. Our first result gives an explicit solution 
of the T-optimal design problem in the case b = 0i jn -i = and — as a by- 
product — proves the conjecture raised in Atkinson (2010). 

Theorem 3.1. A design £ is T-optimal for discriminating between mod- 
els (2.3) and (2.4) with 8\ n -\ = on the interval [—1, 1] if and only if it can 
be represented in the form £ = (1 — a)£i + 0^2 , where a £ [0,1], the mea- 
sures £1 and £2 ar & defined by 

(3.3) & = (*?> - «°>Y ? 2 =(- 4 »<°> - -«°> 

w ■ • ■ oj* r s V u* ■■■ wT 



* 2 ■ 2 ( i7T \ * 2 2 f i7T 



and the weights and support points are given by 
(3.4) 

and (3.2) for b = 0, respectively. 



* 1 

u n = ~ 
n 
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Proof. It was proved by Dette and Titoff (2009) (see Theorem 2.1) 
that any T-optimal design on the interval [—1,1] for discriminating between 
the polynomials X)?=o ®i,j x ^ ano - 

n-2 

r ]1 (x,9 1 ) = ^29 ld x j + 6 ln x n 

3=0 

(note that 0\ n -\ = 0) is supported at the set of the extremal points 

i=(xe[-i,i]|f(i)= sup |^*(i)|}, 
L te[-i,i] J 

where tp*(x) = rji(x,6i) - ^?= 02j2^ and 



j 

n-2 

(3.5) 02 = (^2,o,---,^2,n-2) T = argmin sup 

feel™- 1 ase[-i,i 



rti{x,6 x ) -^02,j 

j=0 



is the parameter corresponding to the best approximation of r/i(x,6i) with 
respect to the sup-norm. By a standard result in approximation theory 
[see Achiezer (1956), Sections 35 and 43] it follows that the solution of 
the problem (3.5) is unique and given by ip*(x) = 9i )n 2~^ n ~ 1 ^T n (x), where 
T n (x) = cos(narccosx) is the nth Chebyshev polynomial of the first kind. 
Note that T n (x) is an even or odd polynomial of degree n with leading coef- 
ficient 2 n ~ 1 [see Szego (1975)]. The corresponding extremal points are given 
by s = *i(0) = -l, Xi = t?(0) = -cosf , i = l,...,n-l, x n = t* n (0) = l. 

Now it follows from Theorem 2.2 in Dette and Titoff (2009) that a de- 
sign is T-optimal if and only if it satisfies the system of linear equations 

(3.6) / 4>* {x)x k dC {x) = 0, k = 0,...,n-2. 

[Note that in the case of linear models the necessary condition in Theo- 
rem 2.2 in Dette and Titoff (2009) is also sufficient.] Therefore for proving 
that £J = £i is a T-optimal design, it is sufficient to verify the identities 

(3.7) f r(x)dei(x)=0i,n2-( n - 1 \-irj2(- 1 ) ix i»i=° 

J i=l 

{k = 0, 1, . . . ,n — 2), which will be done in the Appendix. In a similar way 
we can check that the design £| m (3-3) is a T-optimal design. Note that 

su PP(£i) U SU PP(^2) = = — cos^— i\ i = 0, . . . ,n\ = A, 

because i*_j(0) = —t*(0). Moreover, (3.6) defines a system of linear equa- 
tions of the form Flu = for the vector to = (ojq, . . . ,co n ) T of the T-optimal 
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design £*, where the matrix F is given by F = ((— l) l a^)£z"o '"'^~ 2 £ ]K n_lxn+1 
and has rank n — 1. Additionally, the components of the vector cj sat- 
isfy YH=o w « = 1 • Therefore the set of solutions has dimension 1 . Because 
the vectors of weights corresponding to the designs £* and £?j are given 
by cjW = (0,o;*, . . . ,w*) T and oj^ = (cj* , . . . , u*, 0) T and are therefore lin- 
early independent (note that lo* > 0, i = 1, . . . , n), any vector of weights cor- 
responding to a T-optimal design must be a convex combination of u/ 1 ' 
and u)( 2 \ Consequently, any T-optimal design can be represented in the 
form £ = (1 — + ctQ, which proves the assertion of Theorem 3.1. □ 

Note that the T-optimal design is not unique in the case 6 = 0. On the 
other hand, the T-optimal designs are unique, whenever 9\ n _i ^ 0, and, if 
the ratio |#i,n-i/#i,n| is not too large, the T-optimal designs can also be 
found explicitly as demonstrated in our following result. 



Theorem 3.2. If the parameter b = 9x n —\/6\ n satisfies (3.1), then there 
exists a unique T-optimal design on the interval [—1,1] for discriminating 
between models (2.3) and (2.4). For positive b this design has the form 



(3.8) C 



q(b) ••• t*(b) 



where the points t*(b) and weights w*{b) are defined in (3.2) and (3.4), 
respectively [note that t*(b) > — l,i*(6) = 1]. The T-optimal design for neg- 
ative b has the form 



<(6) ••• -t\(b) 



[note that -t*(6) = -l,-t\{b) < 1/. 



Proof. We consider the case < b < n(l - cos(^))/(l +cos(^)) where 
direct calculations show that the points t*(b), i = 1, . . . ,n, are contained in 
the interval [—1,1]. Moreover, these points are the extremal points of the 
polynomial 



(3.9) 



CnTn 



-x — b/n 
1 + b/n 



b 

1 + - 

n 



where T n is the Chebyshev polynomial of the first kind. For later purposes 
we note that the coefficient of x n_1 in this polynomial is equal to 



(3.10) 



E 

1=1 



b\ b 

1 + - )Ui+- 

n I n 



b, 
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where u\, . . . , u n are the roots of the polynomial T n (x), that is, U{ = cos(^j^-7r) 
{% = 1, . . . , n), Y17=i u i = It can be shown by a standard argument in ap- 
proximation theory [see Achiezer (1956), Sections 35 and 43] that 9i n ip*(x) 
with 

/ -x-b/n 

ip (x) = c n T n 

\ 1 + o/n 

is the unique solution of the extremal problem 

n-2 



mm sup 



where r}x(x,0{) = Y2j=o@i,j xJ ■ Therefore by Theorems 2.1 and 2.2 in Dette 
and Titoff (2009), a T-optimal design is supported at the n extremal points 
ti(b), . . . ,t^(b) [note that we use b < ntan 2 (^-) at this point, which implies 
\t*j (b)\ < 1; j = 1, . . . ,n] and the weights are determined by (3.6). Because the 
set of extremal points is given by A = {t\(b), . . . , i n (&)}, this system reduces 
to 



(3.11) J2t* k (b)(-iyoj* = 0, k = 0,1,..., n-2, 

=1 



and we will prove in the Appendix that the weights given in (3.4) define 
a solution of (3.11). Therefore the design £* specified in (3.8) is a T-optimal 
design for < b < n(l — cos7r/n)/(l + cos7r/n). Since the function ip*(x) is 
unique, any T-optimal design is supported at the points t\(b), . . . ,t^(b) [see 
Theorem 2.1 in Dette and Titoff (2009)]. By Theorem 2.2 in the same refer- 
ence, it follows that the weights of any T-optimal design satisfy the system 
of linear equations (3.11) with uj* = cjj and Y17=i u i = 1- Since ijj*(t*(b)) = 
(— 1)* (i = 1, . . . , n) we can rewrite this system as 

(3.12) Fuo = e n , 

where lo = (u>\, . . . ,io n ) T is the vector of weights, the last row of the ma- 
trix F is given by (1, . . . , 1) and corresponds to the condition X^=i w « = 1) 
e n = (0, . . . , 0, 1) T € W 1 denotes the nth unit vector and the columns of the 
matrix F are given by 

ai = (-i) f (i, t?(6), . . . , (t*(b)) n - 2 ,r(tm)) T , i = 1,2, . . . , n . 

The remaining assertion of Theorem 3.2 follows if we prove that detF ^ 0, 
which implies that the solution of (3.12), and therefore the T-optimal design, 
is unique. For this purpose assume that the opposite holds. In this case the 
rows of the matrix F would be linearly dependent, and there exists a vector 
h = (hi, . . . , h n —\, l) r such that afh = 0,i = 1,2, . . . ,n. But the function 
k(x) = (l,x,. . . ,x n ~ 2 (x)) T h is a polynomial of degree n with coefficient 
of x n_1 given by b. Since aih = k(t*(b)) = this polynomial has roots at the 
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Table 2 

The critical values = ntan 2 (^-) for various values n 6 N 



n 


3 


4 


5 


6 


7 


8 


9 


10 


K 


1 


0.6864 


0.5280 


0.4306 


0.3646 


0.3168 


0.2801 


0.2509 



points t*(b), moreover 

Ew)=-»-eKM$)=-* + i + £. 

i=i i=i v 7 v 7 

However, by (3.10) the sum of the roots must equal —b by Vieta's formula. 
This contradiction proves that detT ^ 0. Therefore the system of equations 
in (3.12) has a unique solution, which means that the T-optimal design is 
unique. 

The case of negative b is considered in a similar way, and the details are 
omitted for the sake of brevity. □ 

The critical values 6* = ntan 2 (^-) for various values of n G N are displayed 
in Table 2. Theorems 3.1 and 3.2 give an explicit solution of the T-optimal 
design problem for discriminating between a polynomial regression of degree 
n — 2 and n, whenever \b\ = |#i,n-i|/|0i,n| < b n . In the opposite case the 
solution is not so transparent and will be discussed in the following section. 

4. T-optimal designs for large values of |6|. In this section we consider 
the case \b\ >ratan 2 (^), for which the T-optimal design cannot be found 
explicitly. Therefore we present a numerical method to determine the op- 
timal designs. The method was described by Dette, Melas and Pepelyshev 
(2004) in the context of determining optimal designs for estimating individ- 
ual coefficients in a polynomial regression model [see also Melas (2006)], and 
for the sake of brevity, we only explain the basic principle. For this purpose 
we rewrite the function fj in (2.5) as 

(4.1) fj(x, a,b) = a + a x x H h a n - 2 x n ~ 2 + 6>i n _i(x n-1 + bx n ), 

where 6 = 1/6 = 9\ n /6\ n —\. Note that for fixed 6, the T-optimal design is 
independent of the parameter 6q n _i and that the choice 

6G --cot 2 f — ),icot 2 f — ) 
n \2n J n \2n J 

corresponds to the case |6| > ntan 2 (^-) considered in this section. In order 
to express the dependence on the parameter 6, we use the notation t*(b) for 
the support points and oo* (6) for the weights of the T-optimal design in this 
section. 

The main idea of the algorithm is a representation of the support 
points t*{b) and corresponding weights oj*(b) in terms of a Taylor series, 
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where the coefficients can be determined explicitly as soon as the design is 
known for a particular point b. The algorithm proceeds in several steps: 

(1) Initialization: In the present situation the point b is given by b = 0, 
which corresponds to the situation of discriminating between a polynomial 
of degree n — 2 and n — 1 . For this case it follows from Dette and Titoff 
(2009) that the T-optimal design coincides with the L>i-optimal design. This 
design has been determined explicitly by Studden (1980) and puts masses 

= rri a * the points ij(0) = cos( ^~^| 7r ) (i = 2, . . . ,n — 1) and masses 

oji(0) = 0J n (0) = 2 (n~i) a ^ P°i n ts ti(0) = — 1 and t n (0) = 1. 

(2) T/ie cinol problem: For the constructions of the Taylor expansion we 
now associate to each vector 



(h, ■ ■ -^n-ij^i) • • • ,w n -i) J 



-Kt 2 < ■■■ <t„_i < 1; 



n-l 

> 0, < 1 



a design with n support points defined by 

Wl W2 ■ ■ ■ w n _i 



0J„ 



As pointed out in the previous discussion, there exists a corresponding ex- 
tremal problem defined by 

(4.2) inf sup \bx n + x^ 1 - f r (x)q\ 



xe[-i,i] 

with a unique solution corresponding to the T-optimal design problem under 
consideration, where we use the notation f T (x) = (l,x, . . . ,x n ~ 2 ). 

(3) The necessary condition: For each vector q in (4.2), define vectors 
r = (q,r) and a quadratic form 



d g = ((?V,6) 



H(e~b)=H(q,T,b)=d^M^ T )d q , 



where M(£ r ) is the information matrix of the design £ r for the regression 
model (4.1). It then follows by similar results as in Dette, Melas and Pe- 
pelyshev (2004) that the design £ T * is a T-optimal design for discriminating 
between the polynomials of degree n and n — 2, and the vector q* is a solu- 
tion of an extremal problem (4.2) if the points G* = (q*,r*) G M™" 1 x IA are 
the unique solution of the system 

d 



H(Q,b) 



0. 



e=e* 



such that the inequality \d^*f(x)\ 2 < d^*M{^ T *)d q * holds for all x € [—1, 1]. 
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(4) Taylor expansion of the optimal solution: The function 
j > .][ t 3n-4 j 

b _> e*(6) = (61(5), • • • , e^ n _ 4 (6)) = (g*(6) r , r*(6f ), 



6*: 



which maps the parameter b~£l = [— ^ cot 2 (^), ^ cot 2 (^)] to the coordi- 
nates of the best approximation q* (b) and the support points t*(b) and 
weights oj*(b) of the T-optimal design, is a real analytical function. The 
coefficients in the corresponding Taylor expansion, 

oo 

e*(b) = e*(b ) + ^2e*(j,bo)(b-b y 

i=i 

in a neighborhood of any point &o G J, can be calculated by the recursive 
formulas 



e*( s + i,6 ) = -7 7 ^^ 1 (&o)(^) S+1 ^ s) (5),fo) 



5=6o 

S = 0,l,2,..., 



where 



%)(b) = %)(bo) + £ Q *(j~ b °)(b ~ boY, 

3=1 

g (e,b) = ^H(e,b), 

J ^)=( 7 ^ V H(@,b)) 

ydQidBj J e=e »(6 ) 

We can use this procedure to calculate the T-optimal design for discrimi- 
nating between polynomials of degree n and n — 2 in the cases which are 
not covered by Theorems 3.1 and 3.2. We illustrate the methodology in the 
following example. 



Example 4.1. Consider the T-optimal design problem for a model of 
degree 5 and a cubic polynomial model. Note that for n = 5, we have 
ntan 2 (^-) ~ 0.528. Therefore if b € [0,0.528], a T-optimal design is given 
by Theorem 3.1, that is, 

f * = (h{b) t 2 (b) t 3 (6) t 4 (b) i\ 

iT ^0.038 0.138 0.262 0.362 ± J ' 
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Fig. 1. The support points (left panel) and weights (right panel) of the T -optimal de- 
sign for discriminating between a polynomial of degree 3 and 5 for various values of 
b = 1/6 € [0,1.894]. 

In order to construct the T-optimal design on the interval [0.528, oo], we 
introduce the notation 6= 1/6 € [0,1.894]. With the results of the previ- 
ous paragraph we obtain a Taylor expansion for the interior support points 
t '^(b), 1 3(b), t%(p) an d weights u;*(&),o;|(5),u;3(5),cj|(&) of the T-optimal de- 
sign for discriminating between a cubic and a polynomial of degree 5 where 
b = 0\n/Q\n~i- By the results of Studden (1980), the vector of support points 
and weights corresponding to the center of the expansion at the point Bq = 
is explicitly known; that is, 

( i *(0),i*(0),^(0),^(0),...,^(0)) = (-^,0,^,|,i,i4). 

At the first step we use a Taylor expansion at the point 60 = to determine 
the T-optimal design for b £ [0,0.4]. When we have found the vector O*(0.4) 
we construct a further Taylor expansion at the point 60 = 0.4, and this 
process is continued in order to determine the vector Q*(b) for any value 
b € [0,1.894]. The support points and weights are depicted in Figure 1 as 
a function of the parameter b = l/b = 0i n /6i n -i. Note that in all cases 6^0 
the T-optimal design for discriminating between a polynomial of degree 5 
and 3 is supported at 5 points. 

5. Concluding remarks and further discussion. In this paper we have 
determined T-optimal designs for discriminating between two rival polyno- 
mial regression models of degree n — 2 and n. To the best of our knowledge 
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these results provide the first analytic solution of a T-optimal discriminat- 
ing design problem with an arbitrary number of parameters in the regression 
model. 

It should be pointed out that the results depend on the ratio of the coef- 
ficients of the terms x n_1 and x n in the polynomial of larger degree, which 
is a well-known feature of the T-optimality criterion. Therefore the designs 
derived here are local in the sense of Chernoff (1953). Usually locally optimal 
designs serve as a benchmark for commonly used designs, as demonstrated 
in the example in the Introduction. Moreover, locally optimal designs form 
the basis for more sophisticated design strategies, which require less knowl- 
edge about the model parameters such as Bayesian or standardized max- 
imin optimality criteria [see Chaloner and Verdinelli (1995) or Dette (1997), 
among others] . This extension was already mentioned in the pioneering work 
of Atkinson and Fedorov (1975a, 1975b) and we conclude this paper with 
a brief discussion of a first explicit result on maximin T-optimal designs for 
the polynomial regression models. 

To be precise, consider the situation where the ratio b = Oi t n-i/Qi,n cannot 
be exactly specified, but prior knowledge suggests that b&I for some interval 
7cR. Without loss of generality, assume 6\^ n = 1; then following Atkinson 
and Fedorov (1975a), a maximin optimal discriminating design maximizes 
the expression 

(5.1) inf inf / ( x n + bx n ~ l + V 9 2 tsA d£(x). 

beie 2 eM.^J_ i y ' 3 ) 

The following result provides a solution of this optimal design problem for 
specific intervals Jcl. 

Theorem 5.1. (a) If I = M, the maximin T-optimal discriminating de- 
sign is given by 



(5.2) e 



AIM 



'o 


t* 


J.* J.* 

" L n-l L n 


1 


1 


1 1 


2n 


n 


n 2n 



where the support points are defined by 



In — % \ 

= cos I 7r , i = 0, . . . , n. 

\ n J 

(b) Assume that I = (— oo, — bo] or I = [bo, oo). If bo > 0, then the maximin 
T-optimal discriminating design coincides with the T-optimal discriminating 
design determined in Sections 3 and 4 for the value b = bo. 

In particular, if bo = 0, then all designs specified in Theorem 3.1 are max- 
imin T-optimal discriminating designs. 
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Proof. In order to prove part (a), note that for I = R criterion (5.1) 
reduces to 

SUP M ! ( X n + hn-l* j ) <%(x), 

which corresponds to the T-optimal discriminating design problem for a poly- 
nomial of degree n and n — 1. By the results in Dette and Titoff (2009), the 
solution of this problem coincides with the Di-optimal design, which is given 
by (5.2) [see Studden (1980)]. 

For a proof of part (b), observe that 

n-2 \ 2 



supinf sup / \x n + bx n 1 + ^ 02 jX j \ d£{x 
H ^fl 2 eR»-W-i\ T^n ' / 



j=0 

1 / n—2 \ 2 

= infsup sup / I x n + bx n - 1 + y"e 2j x j ) d£(x) =: inf R(b), 

where the last equality defines the function R in an obvious manner. We 
now consider the case / = [bo, oo) with bo > and show that the function R 
is increasing on R + , which implies 

(5.3) mfR(b) = R(b ) 

and proves the assertion for the case / = [6o,oo). Recall the definition of 
b* = n tan 2 (-7r/2n) in (3.1); then the proof of Theorem 3.1 shows that for 
6(0,6*] 

6\ 2n 1 



which is obviously increasing with respect to the argument b. If R would 
be not increasing on the remaining region R + \ (0,6*], then there would 
exist real numbers 62 > 61 > 6*, such that R(b±) = R(b 2 ) with corresponding 
extremal polynomials 

L l (x)=x n + b i x n - 1 +qff(x), i = l,2, 

where f(x) = (l,x, . . . ,x n ~ 2 ) T and 

% = argmin f (x n + b iX n ~ l + q T f{x)f d£{x). 

gSR™- 1 J-l 

This yields 

sup \L 1 (x)\= sup |L 2 (x)| = \/R{bi) = ^jR(b 2 ). 

xe[-i,i] xe[-i,i] 

By the discussion in Section 4, the polynomials L\,L 2 can be chosen such 
that they coincide at the boundary points of the interval [—1, 1] (note that 
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for b> b* the support of the optimal discriminating design always contains 
both boundary points —1 and 1). Therefore a simple argument shows that 
there exist n — 2 other points in the interior of the interval (—1, 1), where 
the polynomials must coincide. Consequently, L\(tj) = Lzitj) for n points 
ti,...,t n € [—1,1], which shows that the polynomials are identical. This 
yields b\ = 62, and because of this contradiction the monotonicity of the 
function R has been established, which proves (5.3) and part (b) in the case 
I = [60,00). The remaining case / = (— 00, — 60] can be proved by similar 
arguments, and the details are omitted for the sake of brevity. □ 

Theorem 5.1 provides the solution to maximin T-optimal discriminating 
design problems for specific intervals I C K. In particular, it identifies the 
worst CBjSC clS ct boundary point of the interval under investigation using the 
monotonicity of the criterion with respect to b. This property, which appears 
in many minimax- or maximin optimal design problems, has been criticized 
by Dette (1997). This author recommends Bayesian or standardized max- 
imin optimality criteria, which reflect the different sizes of the optimality 
criteria for different values of b in a more reasonable way. The determina- 
tion of T-optimal discriminating designs with respect to these criteria is 
substantially harder and a challenging problem for future research. 

APPENDIX: PROOF OF IDENTITIES (3.7) AND (3.11) 

Note that the identities in (3.7) and (3.11) can be written in the form 

n 

(A.l) ^i* fc (&)(-l) 4 o;* =0, k = 0, 1, . . . ,n - 2, 

i=i 

where t*{0) = cos(— ) = Xj. We will prove that these equalities hold for any 
real number b. Since 

(A.2) t* k {b) = J2 a j cos (—\ i = 0,l,...,n,/c = 0,l,...,n-2, 

for some coefficients dj = djib) (j = 0,1, ... , k) the identities in (A.l) follow 
from 

(A.3) 2^-1)* cos )w* = 0, fe = 0,l,...,n-2. 

i=i \ n J 

In order to prove (A.3), consider first the case k = 0, n = 2s for some s, where 
the left-hand side of (A.3) reduces to 



i>.*(-D< - 1 g [(1 - »(£)) (-D' + (1 + »(!)) <-!>' 
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+ (-l) s + l 



1 

n 



s-l 



^2(-i)*+(-ir+i 



which proves (A. 3). If k = 0, n = 2s + 1, we get 



£«?(-!)* = - 

i. — < <n 



1=1 



E 



17T 



1-COB - (-1) 



11 



i + C os(|))(-iy 



+ (-i) 



2E«»( r |)(- 1 ) i+1 - 1 



1=1 



cos[vr(l + 2(n + l)s)/(2n)] 



cos(7r/ (2n)) 

l cos((2s + l)7r/2) _ q 
n cos(-7r/(2n)) 

where the third identity follows by standard results for trigonometrical sum- 
mation [see, e.g., Jolley (1961), formula (428)]. This proves (A. 3) for the case 
k = 0,n = 2s + 1. Now consider the case of even n,n = 2s for some odd s, 
s = 21 — 1 and k of the form k = 2(2r — 1). In this case the left-hand side of 
(A. 3) reduces to 

-l 



E 

.i=i 



1 " COS ©) + ( 1 + C ° S (^ 



, kiir\ 



/,'TT 



+ (-l) S COs( ^ ) +C0s(/C7r) 



s-l 



2 E(- 1 ) l -(^)+(- 1 ) 



cos 



(t) 



+ cos(fc7r) 



i 



COS 7T + COS 1 (fc + 2s 

1 4s / V 4s 2 



+ 2 



-{(-l) + (-l) 2 ^ 1 + 2} = 0, 



n 



where we have again used well-known results on trigonometric summation 
[see Jolley (1961), formula (428)]. Therefore we obtain equality (A. 3) in the 
case n = 2s, s = 21 — 1 and k = 2(2r — 1). The other cases can be proved in 
a similar way, and the details are omitted for the sake of brevity. 
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